Studies of Vibrational Properties in Ga Stabilized 5-Pu 
by Extended X-ray Absorption Fine Structure 
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Temperature dependent extended x-ray absorption fine structure (EXAFS) spectra were measured 
for a 3.3 at% Ga stabilized Pu alloy over the range T= 20 - 300 K. EXAFS data were acquired 
at both the Ga K-edge and the Pu Lm-edge. Curve-fits were performed to the first shell interac- 
tions to obtain pair-distance distribution widths, a, as a function of temperature. The temperature 
dependence of a(T) was accurately modeled using a correlated-Debye model for the lattice vibra- 
tional properties, suggesting Debye-like behavior in this material. Using this formalism, we obtain 
pair-specific correlated-Debye temperatures, O c d, of 110.7 ±1.7 K and 202.6 ±3.7 K, for the Pu-Pu 
and Ga-Pu pairs, respectively. The result for the Pu-O c d value compares well with previous vibra- 
tional studies on <5-Pu. In addition, our results represent the first unambiguous determination of 
Ga-specific vibrational properties in PuGa alloys, i.e, c d for the Ga-Pu pair. Because the Debye 
temperature can be related to a measure of the lattice stiffness, these results indicate the Ga-Pu 
bonds are significantly stronger than the Pu-Pu bonds. This effect has important implications for 
lattice stabilization mechanisms in these alloys. 

PACS numbers: 61.10.Ht, 61.66.Dk, 63.20Dj 
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I. INTRODUCTION 

Plutonium in its elemental form presents a complicated 
picture of phase stability!!! Indeed, Pu may adopt one of 
six crystallographically different phases (a , j3, 7 , 6 , 8', 
and e) upon heating from room temperature to its melt- 
ing point of 913 K. This complex phase behavior has been 
explained, in part, by the behavior of the Pu 5/ orbitals 
which fluctuate between itinerant and localized behavior. 
In this heuristic framework, elemental Pu represents the 
transition point along the actinide series from the delo- 
calized electronic nature of the early actinides (Ac-Np) to 
localized, lanthanide-like /-orbital character observed for 
the heavier actinides.cl Theoretical calculations currently 
support the hypothesis that the 5/ orbitals are delocal- 
ized for elements lighter than Pu in the.-actinide series 
and localized for those heavier than PuB This view of 
the relationship between electronic and crystallographic 
structure is reinforced by comparing the complex phase 
behavior and structures observed for pure U, Np, and 
Pu metals as opposed to the relatively simple structures 
observed for the elements Am and beyond. 

The structural and electronic relationships that exist 
between the a and 6 phases have received much attention 
due to some unusual observations. Pure face-centered cu- 
bic (fee) 5-Pu is stable from 593 to 736 K, and exhibits 
a 25% increase in volume relative to that of the ground 
state phase,™ monoclinic a-Pu. However, it is well known 



that the fee structure can be stabilized down to ambient 
temperature by the addition of small amounts (~3-9 at%, 
) of alloying elements such as Al, Ga, In, Sc, and Ce.cl 
However at lower impurity atom concentrations, the S 
phase converts directly to the a form upon coolingj-pos- 
sibly through a martensitic phase transformation.™]! 

Other than a tendency to form trivalent cations, rela- 
tively little is known^about the mechanism of these so- 
called "5-stabilizers" .EH Band structure studies have sug- 
gested that Al, Ga and Sc impurities diffuse the 5/ bands, 
thereby removing their involvement in bonding and lead- 
ing to the stability of-a, more common, (i-bonded transi- 
tion metal-like phasejj Becker et al. have employed an 
ab initio LDA approach to study the lattice relaxation in 
a Pu3iGa supercell cluster, and find evidence of a lattice 
contraction around the Ga site which relax.es the bonding 
constraints fat the neighboring Pu atoms£j Other theo- 
retical worksEi poist to a substantial level of 5/ localiza- 
tion in the <5-phasefc3 and have speculated on a Kondo-like 
model for the involvement of the 5/ electrons£3 

Analysis of the vibrational properties in these materi- 
als is also important for understanding 5-phase stabliza- 
tion. Unfortunately, measurement of the phonon dis- 
persion using conventional inelastic neutron scattering 
techniques is troublesome given the difficulty in growing 
high-quality single crystals. An alternative approach is 
to evaluate the Debye temperature, Qp, using various 
structural techniques. Ledbetter et alXJ employed ultra- 
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sonic wave measurements of elastic constants for 3.3 at% 
Ga 6-Pvl and calculated a Onvalue of 115 K. The thermal 
behavior of 5.0 at% Al 5-Pu was studiaeLby temperature- 
dependent neutron powder diffractioncJ yielding a sim- 
ilarly low value of 0d=132 K. More recently, Lynn et 
al.£3 studied a 3.6 at% Ga <5-Pu sample using neutron- 
resonance Doppler spectroscopy, a technique that can de- 
termine element specific values for 0d • The experiments 
determined a Pu-specific value of 0d=127 K, and as- 
signed a Ga specific value of Od=255 K, although with 
relatively large errors (±22 K). As a test of the Debye 
model, one can check the assumption of equal force con- 
stants for th e Pu and Ga sites by comparing the Ga- 
Od to the y mp n /mG a weighted value of 236 K derived 
from the Pu-0Dvalue. Thus, this comparison does not 
exclude the possibility that the Ga atoms experience a 
stiffer force field compared with the Pu atoms. 

In this article, we present temperature dependent EX- 
AFS (extended x-ray absorption fine-structure) spectro- 
scopic results for 3.3 at% Ga <5-Pu as a means of dis- 
cerning differences in the vibrational character of the Ga 
and Pu sites. Previously, EXAFS studies on Ga stabi- 
lized <5-Pu have focused on isothermal measurements at 
either the Ga K- or the £jjuL«i-edges, and have revealed 
some important effects. In general, these studies 
indicate that the Ga atoms reside in their expected fee 
lattice positions although there is an appreciable lattice 
contraction observed for the Ga-Pu bonds (~3-4%). Sur- 
prisingly, the contraction in Ga-Pu bonds is significantly 
larger than the collapse calculated theoretically!!^ or the 
contraction expected from Vegard's law based on a sim- 
ple substitutional alloy. In addition, EXAFS data show 
increased disorder for the Pu-Pu near neighbor interac- 
tions. These results indicate that there are unexpected, 
site-specific lattice effects occurring in these materials. 

The outline of the paper is as follows: details of sam- 
ple preparation and EXAFS experimental setup and data 
analysis are discussed in Sec. II. The results of curve- 
fitting analysis and modeling the cr(T) temperature de- 
pendence using the correlated Debye model are presented 
in Sec. III. A discussion of the results and their relation 
to other 0d studies is presented in Sec. IV, and the 
conclusions are given in Sec. V. 



II. EXPERIMENTAL DETAILS 

A. Sample Preparation 

A ^6 fim thick 239 Pu foil (3.3 at % Ga) was prepared 
from 20 year old material by melting to remove accumu- 
lated helium, followed by subsequent annealing, cutting, 
and rolling. The foil was further homogenized at 450 °C 
for ^100 hrs to ensure that single-phase, S-Pu was pro- 
duced. Transmission x-ray diffraction was performed at 
LLNL and confirmed the presence of the fee phase, with 
no significant amounts of other Pu phases present. In 
preparation for EXAFS analysis, the foil was electropol- 



ished to remove any accumulated oxide material on the 
surface. The sample was then encapsulated under ar- 
gon using a specially designed, triple containment x-ray 
compatible cell manufactured by Boyd Tech. The first 
level of containment consisted of coating the sample in 
a thin film of liquid polyimidc solution that was allowed 
to air dry directly on the sample. The foil was then 
mounted onto an aluminum frame and sealed within two 
additional layers of x-ray transparent Kapton windows 
(0.010" thick). The windows were clamp-mounted onto 
the aluminum body with stainless steel window frames 
and using indium wire as a vacuum seal material. The 
triple-contained sample was subsequently mounted in an 
open cycle liquid helium flow-cryostat for variable tem- 
perature EXAFS measurements. Temperature measure- 
ment errors are within ~1 K, and are stable within ^0.2 
K. 



B. EXAFS Data Acquisition and Analysis 

Plutonium Lm- and gallium if -edge x-ray absorption 
spectra were collected at the Stanford Synchrotron Radi- 
ation Laboratory (SSRL) on wiggler beamline 11-2 under 
normal ring operating conditions using a nitrogen-cooled 
Si (220), half-tuned, double-crystal monochromator oper- 
ating in unfocussed mode. The vertical slit height inside 
the x-ray hutch was 0.3 mm which reduces the effects 
of beam instabilities and monochromator glitches while 
providing ample photon flux. The Pu Lm-edge spectra 
were measured in transmission mode using Ar-filled ion- 
ization chambers. The Ga K-edge spectra were measured 
in fluorescence mode using a 30-element Ge array solid 
state detector developed by Canberra Industries. The 
detector was operated at ~75 kHz per channel, and the 
signals were digitally processed using the DXP 4C/4T 
developed by X-ray Instrumentation Associates. 

XAFS raw data treatment, including calibration, nor- 
malization, and subsequent processing of the EXAFS and 
XANES (x-ray absorption near-edge structure) spectral 
regions was performed by standard methods reviewed 
elsewhereaO using the EXAFSPAK suite of programs 
developed by G. George of SSRL. Typically, three XAFS 
scans (transmission or fluorescence) were collected from 
each sample at each temperature and the results were av- 
eraged. The spectra were energy calibrated by simulta- 
neously measuring the absorption spectrum for the refer- 
ence samples PuC-2 or Ga2 03. The energies of the first in- 
flection points for the reference sample absorption edges, 
E r , were defined at 18053.1 eV (Pu L m ) and 10368.2 
eV (Ga K). The EXAFS threshold energies, E , were 
defined as 18070 eV and 10385 eV for the Pu and Ga 
edges, respectively. Nonlinear least-squares curve-fitting 
was performed on the fc 3 -weighted EXAFS data using the 
EXAFSPAK programs. 

The EXAFS data were fit using theoretical phase 
and amplitude functions— calculated from the program 
FEFF8.1 of Rehr et al&B All of the Pu-Pu interactions 
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were modeled using single scattering (SS) paths derived 
froiM the model compound, unalloyed <5-Pu, af cc =4.6371 
AM The Ga-Pu SS interactions were modeled by using 
the same model compound structure and replacing the 
central absorbing atom with Ga. An initial series of fits 
was done on the raw Pu Lin and Ga if-edge fc 3 -weighted 
data sets using the expected fee near-neighbor interac- 
tions at 3.28, 4.64, 5.68, and 6.56 A and fixing the coor- 
dination numbers at 12, 6, 24, 12, and 24, respectively. 
The results (especially at low T) confirmed the presence 
of the fee structure and showed no evidence for other un- 
usual structural effects. That is, we observed no phase 
changes (i.e., S -±-&Lor previously postulated metastable 
impurity phases £§Ea This first level of analysis was also 
used to establish values for Sq and AEq by fixing coor- 
dination numbers, N, and allowing Sq, AEq, a 2 , and R 
to vary. The final values used were taken from averag- 
ing over the range of temperatures studied: for Pu Sq = 
0.55, AE = -12 eV; and for Ga, Si = 0.85, AE Q = -10 
eV. 

As a result of the preliminary analyses, all of the sub- 
sequent fits described here focussed on isolating the be- 
havior of the first shell Pu-Pu and Ga-Pu interactions 
in a highly constrained manner. Thus fits were done on 
Fourier-filtered data using the same fixed S^ and AE 
values for all temperatures, along with a fixed coordi- 
nation number of iV=12 for the first shell. The ability 
to fix AEq, N, and 5q helps to avoid correlation prob- 
lems between the fit parameters and to establish more 
consistently any changes in a 2 and R that may occur as 
functions of temperature. 



III. RESULTS 

A. EXAFS Raw Data and Curve-Fitting 

Figures 1A and IB show the raw fc 3 -weighted Pu EX- 
AFS L\i\ data and the corresponding Fourier transforms 
(FT) for the 5-Pu sample as a function of temperature. 
The FT represents a pseudo-radial distribution function 
and the peaks are shifted to lower R values compared to 
real interaction values as a result of the phase shifts asso- 
ciated with the absorber-scatterer interactions (^0.1-0.2 
A for Pu-Pu). As the sample is cooled from 293 K to 20 
K, the EXAFS scattering amplitude increases systemat- 
ically due to decreased thermal motion of the atoms in 
the lattice. As a result, the effective measurable fc-range 
increases substantially at lower temperature. This ef- 
fect is equally visible in the corresponding FTs. As the 
temperature is lowered, the intensity of the FT peaks 
increases dramatically, and the spectra reveal a pattern 
consistent with that expected for a fee lattice. The first 
shell Pu-Pu peak seen at ~3.1 A corresponds to 12 Pu 
near neighbors in the fee structure and increases in max- 
imum peak height while also narrowing with decreasing 
temperature. The second, third, and fourth shell peaks 
at ~4.4, 5.4 and 6.3 A which correspond to real inter- 
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FIG. 1: Pu L m -edge fc 3 -weighted EXAFS data (A) and the 
corresponding Fourier transforms (B) for the S-Pu sample as 
a function temperature. Data were acquired in transmission 
mode and Fourier transformed over the ranges shown in the 
figure. 



actions at 4.64, 5.68, and 6.56 A are clearly affected by 
thermal effects and become distinguishable only below 80 
K. 

The Ga if -edge raw fc 3 -weighted EXAFS and corre- 
sponding FTs are shown in Figures 2 A and 2B. This 
series of spectra qualitatively exhibit the same thermal 
behavior observed in the Pu Lm-edge EXAFS data, that 
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FIG. 2: Ga iC-edge fc 3 -weighted EXAFS data (A) and the 
corresponding Fourier transforms (B) for the (5-Pu sample as 
a function temperature. Data were acquired in fluorescence 
mode and Fourier transformed over the ranges shown in the 
figure. 



is, increased scattering amplitude along with detection 
of more distant neighbors are apparent with decreasing 
temperature. The spectra are dominated by the peak 
at ~3.0 A, attributed to the first shell Ga-Pu interac- 
tions (7V=12). At lower temperatures, the FTs reveal 
second, third, and fourth shell Pu neighbor interactions 
at ^4.4, 5.4, 6.3 A, respectively. In contrast to the be- 



TABLE I: Pu L m and Ga if-edge EXAFS single shell curve 
fitting results. The first shell interactions were isolated by 
Fourier transforming over the data fc-ranges displayed in Fig- 
ures 1A and 2A, and back-transforming over the range R=2-4 
A. Several parameters were held fixed as follows: for Pu and 
Ga, AT=12; for Pu, Sq = 0.55 and AE = -12 eV; and for Ga, 
S 2 = 0.85 and AE = -10 eV. 



Sample 


Pu-Pu 


shell" 


Ga-Pu shell 


Temp (K) 


fl(A) 


a 2 (A 2 ) 


R (A) 


a 2 (A 2 ) 


20 


3.290 


0.00280 


3.160 


0.00275 


50 


3.292 


0.00381 


3.160 


0.00332 


80 


3.293 


0.00511 


3.160 


0.00395 


110 


3.301 


0.00680 


3.160 


0.00512 


140 


3.301 


0.00905 


3.160 


0.00609 


170 


3.306 


0.01039 


3.163 


0.00764 


210 


3.316 


0.01275 


3.165 


0.00953 


293 


3.297 


0.01821 


3.171 


0.01264 



"Errors in R and a 2 are estimated to be ±0.005 A and ±10% 
based on EXAFS fits to known model compounds, cf. Ref. |2l| 



havior observed in the Pu EXAFS, the detection of these 
longer range interactions is retained up to higher temper- 
atures (ca. 110 K). At this level of analysis, this obser- 
vation suggests that the Ga-Pu interactions experience 
less thermal broadening than the corresponding Pu-Pu 
interactions, which is especially surprising given the de- 
creased reduced mass of the Ga-Pu pair relative to the 
Pu-Pu pair. 

Following initial inspection of the data (see Section II) 
we chose to focus our analyses on the isolated first shell 
contributions, in part, due to the high signal-to-noise 
obtained relative to the more distant shell interactions. 
Thus by curve-fitting the Fourier-filtered first shell com- 
ponents from the Pu Lm and Ga K EXAFS, the local 
vibrational temperature dependence may be investigated. 
The first shell interactions were isolated by Fourier trans- 
forming over the data /c-ranges shown in Figures 1A and 
2A, and back-transforming over the range R=2-i A. 

The curve-fitting results summarized in Table Q serve to 
identify the different temperature dependent and static 
structural effects for the Pu and Ga sites in this mate- 
rial. Figure 3 shows the resulting curve fits to the Fourier 
filtered Pu and Ga EXAFS data measured at 20 K in 
fc-space and i?-space. The strong coincidence between 
the data and the fits serves to illustrate the highly re- 
fined nature of the FEFF8.1 calculations as well as the 
appropriateness of the fitting procedure employed. Be- 
cause these fits were highly constrained, any changes in 
the data are assigned to changes in R and a 2 . The first 
effect revealed by these data is the overall contraction 
of Pu atoms around the Ga sites relative to the environ- 
ment around the Pu sites. This "collapse" of about 4% is 
compjaxable-to that observed by earlier EXAFS studies on 
<5-PuOHj'E£I The first shell bond lengths for each element 
appear to remain constant over the measured tempera- 
ture range, at least within the quoted experimental error 
and under the constraints outlined above. 



5 




1 1 II 1 II 

; ( 


1 1 1 1 1 1 

1 * 






Ga 

■ Ar 




V 

Pu 




(x2) 









2 2.5 3 3.5 4 
R(A) 



FIG. 3: Representative single shell curve fits (dotted lines) to 
the Pu Liii and Ga K-edge data (solid lines) measured at 20K 
in (A) fc-space and (B) ii-space. Data were Fourier-filtered 
from the raw spectra over the k and R ranges shown in the 
plots using Gaussian window functions with half-widths of 0.5 
A -1 and 0.05 A, respectively. Pu data shown are multiplied 
by a factor of 2 for the purpose of comparison to the Ga data. 



B. Vibrational Analysis 

The most dramatic effect depicted in this data set is 
the difference in the temperature dependence between 
the Pu-Pu and Ga-Pu Debye- Waller factors, a 2 . The 
values for a 1 increase with temperature consistent with 
greater thermal disorder, as expected. However, accord- 
ing to Table Q, the Pu-Pu first shell shows much greater 
disorder at higher temperature than the corresponding 
Ga-Pu shell. Moreover, this occurs in spite of the fact 
that the Ga-Pu pair has a lower reduced mass. To study 
this effect more carefully, the temperature dependence 
of the Debye- Waller factors was modeled by employ- 
ing the correlMed-Debye Model to determine the Debye 
Temperature. E3 



s (T) = <J 2 taUc + F(T,8, 



cD) 



(1) 



The temperature-dependent part of the Debye- Waller 
factor F(T, 9 C £>) is given within the correlated-Debye 
model by 



ft f 

F{T,9 cD ) = —J P] (uj)coth 



2k B T 
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Ga-8 n 
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202.6 K 
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FIG. 4: Temperature dependence of EXAFS Debye-Waller 
factors for first shell Pu-Pu (upper data) and Ga-Pu (lower 
data) interactions plotted along with fits generated using the 
correlated Debye model. 



Rj is 



3cj 2 
Pi = -J- 



1 - 



sin(ujRj/c) 
toRj/c 



(2) 



where fi is the reduced mass, 9 c d is the correlated Debye 
temperature, and the phonon density of states at position 



in which ujd is the usual Debye frequency and c = uin/kB 
where ks is Boltzmann's constant. The expression in 
brackets of Eq. || takes into account the correlated mo- 
tion of the atom pairs. 

Using this formalism, we obtain the fits shown in Fig- 
ure H and determine pair-specific correlated-Debye Tem- 
peratures, 9 cD , of 110.7 ±1.7 K and 202.6 ±3.7 K, for 
the Pu and Ga sites, respectively. Estimates of the <J 2 tatic 
static parameter are all consistent with no static disor- 
der, and fall in the range of 0.0003 ±0.0003 A 2 . The 
data appear to fit this model quite accurately, indicating 
that the local pair vibrations can be described using the 
correlated-Debye model. However, the higher 9 c d of the 
Ga-Pu pair with respect to the Pu-Pu pair may indicate 
the presence of a local lattice anomaly from a vibrational 
standpoint, which could be related to the observed lattice 
contraction measured around the Ga sites. 



IV. DISCUSSION 

It is instructive to compare the EXAFS 9 c d results de- 
scribed here with those obtained from previous studies. 
Table || summarizes the 9d values for (S-Pu determined 
from our study along with those from earlier works. We 
also make special note of the differences in the methods 
used to determine these values, which is important for 
a proper interpretation. The elastic values of Od are 
obtained by ultasound measurements of "bulk" ..lattice 
wave properties extrapolated down to T=0 K.oESEj The 
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TABLE II: Comparison of Pu and Ga Specific Debye Tem- 
peratures Obtained for <5-Pu by Various Techniques 



<5-Pu alloy 


Pu-9 D 


Ga-6 D 


Method 


3.3 at% Ga 


110.7 




EXAFS 






202.6 


EXAFS 




127 




Neut. Abs." 






255 


Neut. Abs. 




105 




Bulk Mod. 6 




115 




Elastic c 


6.6 at% Ga 


127 




Elastic 1 * 


5.0 at% Ga 


132 




Elastic e 




132 




Neut. Diff/ 



"Neutron Resonance Doppler spe ctroscopy. Ref. [if]. 
'Using the relation from Ref. blland the experimental bulk mod- 
ulus noted in the text. 







'Ultrasonic measurement of 3 principal elastic constants, Cij, Ref. 



Derived from Young's modulus, E, and torsional modulus, G, 
Ref. §. 

e Also derived from E and G. Ref. BCj. 
^Neutron diffraction, Ref. |l5[ 



"Bulk Modulus" value uses a theoretical relationship!! 
where 8d= 41.63^/roB/M to derive Od from an ex- 
perimental value for B=A50 kbar. The Od value from 
neutron diffraction also represents a bulk measurement 
of lattice thermal behavior£j However, since it is de- 
rived from thermal displacement factors specifically at- 
tributed to the Pu atoms, it is referenced as ©dw- The 
©d values pobtained from EXAFS and neutron resonant 
absorptionllj are distinct from those derived from the 
other methods in that they are element specific, yet there 
are important differences to be noted. EXAFS inherently 
measures only phonon modes that involve specific atom 
pairs (i.e., we must consider the reduced mass for the 
pair). In contrast, neutron resonance absorption mea- 
sures the average of all modes associated with a single 
atomic species (i.e., only the mass of the absorbing atom 
must be considered). 

In general the Pu specific values of Od obtained from 
all the techniques are similar in that they are all quite 
low, ranging from 104 to 132 K. It is not our intent 
to discuss these differences since they come from such 
a wide variety of techniques and variability in samples. 
The Ga-spccific Debye temperatures do however require 
more careful discussion, in relation to the Pu-specific val- 
ues. As mentioned in Sec. I, neutron resonant absorp- 
tion found a Ga specific Od value of 255 K ±22 K which, 
compared to the y/ mp u /mG a weighted value of 236 K, 
does not exclude a slightly stiffer force field than the one 
around the Pu sites. Analogously, EXAFS yields a Ga 
specific © c d of 202.6 ±3.7 K which m ay be com pared to 
the expected value of 164 K (i.e., the \/ /j,p u /fJ-Gn weighted 
value). The observation of a significantly higher O c d for 
the Ga-Pu pair relative to that expected for a simple sub- 
stitutional lattice model, coupled with the anomalously 
large lattice contraction around the Ga sites clearly sug- 
gest that the Ga sites reside in a significantly stronger 



force field. The shorter Ga-Pu bonds, in fact, inherently 
imply a strengthening of the Ga-Pu bonds relative to the 
Pu-Pu bonds. 

Regarding factors that contribute to (5-phase stabiliza- 
tion, these results are consistent with the Ga atoms hav- 
ing a significant influence on the local electronic struc- 
ture. One obvious possibility is that charge transfer ex- 
ists between the Ga pd states and the corresponding Pu 
df states. This orbital overlap will directly dictate the 
extent of lattice vibrational distortions and local geome- 
try within the crystal structure, that is, with a preference 
towards a fee lattice in the presence of Ga. Therefore, the 
results of this EXAFS study give the first unambiguous 
evidence for the impact of the vibrational properties on 
the local lattice effects. 

From the results of this study, we may also estimate 
values for "pair-specific bulk moduli." Indeed, accord- 
ing to the empirical relation between bulk modulus, and 
Debye temperature established by Moruzzi et a/.Ej, we 
deduce from our calculated Debye temperatures for Pu 
and Ga values of 498 and 1666 kbar, respectively, for 
the bulk modulus. Since the Pu-Ga alloy we have been 
studying has a composition in the dilute limit, this im- 
plies that Pu is mostly surrounded by Pu atoms whereas 
each Ga site can be viewed as embedded in a Pu matrix. 
In addition, the balk modulus of pure Ga metal as deter- 
mined by theory,E2l 606 kbar, or experiment,^ 613 kbar, 
is quite low in accordance with the low melting point of 
308 K. Hence the dramatic difference between the two es- 
timated bulk moduli provides an equivalent description 
for the increase in strength gained in going from Pu-Pu 
(or Ga-Ga bonds) to Ga-Pu bonds. This unusually large 
increase in stability of a Pu (Ga) matrix by addition of 
Ga (Pu) is reflected in the existence of very stable Pu-Ga 
compounds (i.e., large heats of formation), some of them 
exhibiting congruent melting. It would be interesting to 
extend this work to other alloy compositions to quan- 
tify more appropriately these findings and have a better 
understanding of the unusual synergistic effects due to 
alloying in Pu-Ga and related systems. 



V. CONCLUSION 

New pair specific Debye temperatures for the Ga- 
Pu and Pu-Pu pairs in S-Pu were determined using a 
correlated-Debye model fit to temperature dependent 
EXAFS Debye-Waller parameters. Using this formal- 
ism, we obtain pair specific correlated-Debye tempera- 
tures, O cD , of 110.7 ±1.7 K and 202.6 ±3.7 K, for the 
Pu-Pu and Ga-Pu pairs, respectively. The results for 
the © c d Pu-Pu pair compare well with previous vibra- 
tional studies on <5-Pu. In addition, our results represent 
the first unambiguous determination of Ga-specific vi- 
brational properties, i.e, Ga-Pu c d , in PuGa alloys. 
Because the Debye temperature can be related to a mea- 
sure of the lattice stiffness, these results indicate that the 
Ga-Pu bonds experience a stronger force field than the 
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corresponding Pu-Pu bonds. This effect has important 
implications for lattice stabilization mechanisms in these 
alloys. If the vibrational properties in <5-Pu are dependent 
on the type of impurity, further studies of Debye tem- 
peratures with other "stabilizers" will be important for 
discerning key aspects of the stabilization mechanism. 

Acknowledgments 

This work was performed under the auspices of the 
U.S. Department of Energy (DOE) by the University of 



California Lawrence Livermore National Laboratory un- 
der contract No. W-7405-Eng-48. This work was par- 
tially supported (C. H. B.) by the Office of Basic Energy 
Sciences, Chemical Sciences Division of the U. S. DOE, 
Contract No. DE-AC03-76SF00098. The authors also 
wish to thank John Rehr and Alex Ankudinov for help- 
ful discussions. This work was done (partially) at SSRL, 
which is operated by the Department of Energy, Division 
of Chemical Sciences. 



* Electronic address: allen42@llnl.gov 

1 W. N. Miner and F. W. Schonfeld, in Plutonium Handbook, 
edited by O. J. Wick (The American Nuclear Society, La 
Grange Park, IL, 1990), p. 33. 

2 J. L. Smith and E. A. Kmetko, J. Less-Common Met. 90, 
83 (1983). 

3 J. M. Willis and O. Eriksson, Phs. Rev. B 45, 13879 (1992). 

4 P. Chiotti, V. V. Akhachinskij, I. Ansara, and M. H. Rand, 
The Chemical Thermodynamics of Actinide Elements and 
Compounds (The American Nuclear Society, International 
Atomic Energy Agency, Vienna, 1981), vol. 5, p. 231. 

5 P. H. Adler, G. B. Olson, M. F. Stevens, and G. F. Galle- 
gos, Acta Metall. Mater. 40, 1073 (1992). 

6 T. G. Zocco, M. F. Stevens, P. H. Adler, R. I. Sheldon, 
and G. B. Olson, Acta Metall. Mater. 38, 2275 (1990). 

7 J. R. K. Gschneidener, R. O. Elliott, and V. O. Struebing, 
in Plutonium 1960 (Cleaver-Hume, London, 1961), p. 99. 

8 P. H. Adler, Metall. Trans. A 22, 2237 (1991). 

9 P. Weinberger, A. M. Boring, and J. Smith, Phys. Rev. B 
31, 1964 (1985). 

10 J. D. Becker, J. M. Wills, L. Cox, and B. R. Cooper, Phys. 
Rev. B 58, 5143 (1998). 

11 P. E. A. Turchi, A. Gonis, N. Kioussis, D. L. Price, and 
B. R. Cooper, in Proceedings of the International Work- 
shop on Electron Correlations and Materials Properties, 
edited by A. Gonis, N. Kioussis, and M. Ciftan (Kluver 
Academic and Plenum Publishers, New York and London, 
1999), pp. 531-537. 

12 M. Penicaud, J. Phys.: Condens. Matter 9, 6341 (1997). 

13 S. Meot-Reymond and J. M. Fournier, J. Alloys and Comp. 
232, 119 (1996). 

14 H. M. Ledbetter and R. M. Moment, Acta Metall. 24, 891 
(1976). 

15 A. Lawson, J. Goldstone, B. Cort, R. Sheldon, and E. Foly- 
ton, J. Alloys and Comp. 213, 426 (1994). 

16 J. Lynn, G. Kwei, W. J. Trela, V. W. Yuan, B. Cort, R. J. 
Martinez, and F. Vigil, Phys. Rev. B 58, 11408 (1998). 

17 L. Cox, R. Martinez, J. H. Nickel, S. Conradson, and 
P. Allen, Phys. Rev. B 51, 751 (1995). 



18 P. Faure, B. Deslandes, D. Bazin, C. Tailland, R. Doukhan, 
J. M. Fournier, and A. Falanga, J. Alloys and Comp. 244, 
131 (1996). 

19 N. Richard, P. Faure, P. Rofidal, J. L. Truffier, and 
D. Bazin, J. Alloys and Comp. 271-273, 879 (1998). 

20 T. M. Hayes and J. B. Boyce, in Solid State Physics, edited 
by H. Ehrenreich, F. Seitz, and D. Turnbull (Academic, 
New York, 1982), vol. 37, p. 173. 

21 G. G. Li, F. Bridges, and C. H. Booth, Phys. Rev. B 52, 
6332 (1995). 

22 J. J. Rehr, J. Mustre de Leon, S. I. Zabinsky, and R. C. 
Albers, Phys. Rev. B 44, 4146 (1991). 

23 J. J. Rehr, J. Mustre de Leon, S. I. Zabinsky, and R. C. 
Albers, J. Am. Chem. Soc. 113, 5135 (1991). 

24 F. H. Ellinger, J. Metals 8, 1256 (1956). 

25 S. D. Conradson, Appl. Spectrosc. 52, 252A (1998). 

26 S. D. Conradson, in Los Alamos Science, No 26, Vol. II 
(Los Alamos National Laboratory, Report No. LA-UR-00- 
4100, 2000), pp. 356-363. 

27 E. D. Crozier, J. J. Rehr, and R. Ingalls, in X-Ray Ab- 
sorption: Principles, Applications, Techniques of EXAFS, 
SEXAFS, XANES, edited by D. Konigsberger and R. Prins 
(Wiley, New York, 1988), p. 373. 

28 G. B. Beni and P. M. Platzman, Phys. Rev. B 14, 1514 
(1976). 

29 J. C. Taylor, P. F. T. Linford, and D. J. Dean, J. Inst. 
Metals 96, 178 (1968). 

30 M. Rosen, G. Frez, and S. Shtrikman, J. Phys. Chem. 
Solids 30, 1063 (1969). 

31 V. L. Moruzzi, J. F. Janak, and K. Schwarz, Phys. Rev. B 
37, 790 (1988). 

32 P. E. A. Turchi, A. Gonis, and N. Kioussis, calculated 
for hypothetical fec-based Ga using the FP-LMTO (Full- 
Potential Linear Muffin- Tin Orbital) method, unpublished 
work. 

33 R. W. G. Wyckoff, Crystal Structures, 2nd ed. (Wiley, New 
York, 1962), vol. 1, p. 22. 



